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Abstract 

The QCD inspired potential model persued by us earlier has been 
recently modified to incorporate an additional factor V in the linear 
cum Coulomb potential. While it felicitates the inclusion of standard 
confinement parameter b = 0.183GeV 2 , unlike in previous work, it 
still falls short of explaining the Isgur-Wise function for the B mesons 
without adhoc adjustment of the strong coupling constant . 
In this work, we determine the factor ' d from the experimental values 
of decay constants and masses and show that the reality constraint 
on V yeilds bounds on the strong coupling constant as well on slope 
and curvature of Isgur-Wise function allowing more flexibility to the 
model. 
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1 Introduction 

In recent years, considerable experimental and theoretical efforts have been 
undertaken to understand the physics of hadrons containing a heavy quark 
PQ. The Isgur-Wise function [2] is an important quantity in this area of hadron 
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physics.lt is in this spirit, that this function has been studied in various quark 
models [3] — [12] besides QCD sum rule approach [T3], the MIT bag model [H] 
and the Skryme model p~5| . 

Since one of the basic ingredients of the IW function is the hadron wave- 
function involving heavy quark [3] [12], it is therefore meaningful to test any 
specific QCD inspired quark model by calculating the IW function and study 
it phenomenologically. 

Sometimes back, a specific QCD inspired quark model was proposed by us 
[16J which had later been used to calculate IW function as well [T7J [18j [19]. 
One of drawback of the model is that significant confinement effect could 
not be accomodated in the model[TSJ [T7J [TS] due to perturbative constraints 
coming from using the Dalgarno's method [20]. Only recently [32], standard 
confinement effect b = 0.183Ge^ 2 [21] was accomodated in the improved 
version of QCD inspired quark model brought through the introduction of 
parameter ' d in the potential : V = ~^r- + br + c taking c ~ lGeV as its nat- 
ural scale fixing and Aq = 1 where Aq is an undetermined factor appearing in 
the series solution of the Schrodinger equation [ Eq.(8) of Ref.l9].In earlier 
works [T6l [T71 [18], the unknown coefficient cAq occured in the wavefunction 
was set to zero. 

One of the drawback of work [19]was the adhoc enhancement of strong cou- 
pling constant was needed to take into account of the slope and curvature of 
B,B S and B c mesons. 

In this work, we take an alternative strategy to remove this adhoc en- 
hancement. We use the wavefunction at the origin involving the unknown 
coefficient cA Q and fix it from the experimental values of masses and decay 
constants directly. The reality constraint on cA will then yeild lower bounds 
on the strong coupling constant a s , which would lead to the upper bounds 
on the slope and curvature of the IW function. 

The rest of the paper is organised as follows : Section 2 contains the theory 
of the improved QCD inspired quark model, Section 3 encloses the results 
and in Section 4 we draw conclusion and remarks. 
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2 Theory 

2.1 The Wavefunction 

The spin independent Fermi-Breit Hamiltonian for ground state (I = 0), neglecting 
the contact term proportional to 5 3 is [16j [17] : 

H = H + H' 

V 2 Aa s 

= h or 

2fi 3r 

where a s is the running coupling constant , b is the confinememt parameter 
and c is another parameter whose significance will be cleared later. 
In this work,the a s values are taken from the V -scheme [T8[ [26], [27] as done 
in [19] which are large as compared to those of MS'-scheme.It is necessary as 
large a s values lead to better results for slope and curvature of Isgur-Wise 
function in earlier work [T71 HEJ 119] - 

In earlier work [T7J [18] , large confinement was found to be inconvenient in 
the calculation of slope and curvature of Isgur -Wise function which was a 
limitation of the model. So, our aim has been focussed in the inclusion of 
larger confinement and hence we retain the same choice of b = 0.183GeV 2 
[T9l 121] to investigate whether this approach leads to better results or not. It 
is worth notable that inclusion of c in the analysis [19] allows larger b useful 
, so we don't think another choice of b. 

With ,H = — y~ — p< as the parent Hamiltonian and H' = br + c as the 
perturbed Hamiltonian , we obtain a ground state wavefunction upto the 
first order correction using the Dalgarno method [20J of stationary state per- 
turbation theory as : 



(i) 



VW (r) = N U 4- -f= - e-4 (2) 

where Aq is the unknown coefficient appearing in the series solution of the 
Dalgarno method 

Including the relativistic effect [221 E3] , the wavefunction is : 
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^conf+rel (r) = N' I cA + -J= - ^2^1 (- ) • " (3) 



7ran \ na; 



a 



Here a is given by: 



a 



4/ia s 



(4) 



and 



^-l/l-f (5) 



A" and N' are the normalization constants given by : 



AT 2 = (6) 



1 + _ 3/U&a 3 + ^2^2 + 

and 



'Tin 



r>7-2e 

/2 



AT = = (7) 

r(3-2e)X 1 {) 

where X x is given in APPENDIX A. 

We note that the equations (2),(3),(6)and (7) are obtained from Eq.(4),(6),(5) 
and (7) of Ref. [19] exhibiting explicit dependence of cA in them. 

2.2 Fixing of the coefficient cAq 

The wavefunction at the origin (WFO), is related to the decay constant f p 
and the mass of the pseudoscalar meson M p through the relation 
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Again from equation (2), we have : 



f 2 M 

hMo) r = ^# (8) 
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IV(o)i 2 = iv 2 [ c 2 ^ + -3 + 



2cAr 



nan 



TTCLq 



Using equation(6), we arrive at the quadratic equation for cA : 

A' (cA ) 2 + B' (cA ) + C = 

where 



B' 



A' = 7rag|^ (0) | 2 - 1 
\ip (0) \ 2 - 3{iba 3 ^4\^ (0) 



and 



C=\if>(0) 



1 + 



45/i 2 &X 



(9) 

(10) 

(11) 
(12) 

(13) 



Using the experimental values of f p and M p [25] , we determine \ip (0) | 2 from 
equation(9) which in turn will yeild two solutions for cA in equation (10): 



cAr 



-B' ± VB' 2 - AA'C 
2A' 



(14) 



which will depend on fi,M P ,f P and a s .The solution corresponding to the 
+ve(-ve) sign of equation(15)will be termed as +ve(-ve) solution hereafter.lt 
will be shown numerically that for a given /x, Mp, and fp ,a s reaches the 
minimum value when the following condition is satisfied : 



B' - AA'C 







(15) 



The formalism involving Eq.(5)-(16) is strictly valid only without relativistic 
effect as the wavefunction at the origin with such effect [Eq.(3)] is not well 
defined due to its singularity at the origin. For a subsequent analysis ,we as- 
sume that cA does not deviate significantly from its non-relativistic value so 
that it can be used to calculate the slope and curvature of the IW function 
even without relativistic effect. 
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2.3 Charge radius (slope)and convexity parameter (cur- 
vature) of I-W function 

The Isgur-Wise function is written as [21 E] : 

= l- p 2 ( y -l) + C(y-l) 2 + ... (16) 

where 

V = Vp.v'p (17) 

and and v' being the four velocity of the heavy meson before and after 
the decay. The quantity p 2 is the slope of I-W function at y = 1 and known 
as charge radius : 



The second order derivative is the curvature of the I-W function known as 
convexity parameter : 



d 2 y \y=i 



(19) 



For the heavy-light flavor mesons the I-W function can also be written as 

isiini : 



where 



Anr 2 \ip (r)| 2 cos prdr 



p 2 = 2fi{y-l) 



(20) 



(21) 



Equation (21) holds good for both relativistic and nonrelativistic case. The 
wavefunction ip (r) takes different form for both the cases. Without relativis- 
tic effect, it is given by equation(2) and with relativistic effect it is given by 
(3). 

With the wavefunction(2)in equation(lO) i.e. including confinement only 
the charge radius pl on f and convexity parameter C con f are respectively given 
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by: 

2 fi 2 [24irc 2 A 2 a 5 + 24a 2 + 630p 2 b 2 a 8 + 48cA \f^aJ - 180cA fib\f^ - I80pba 5 ] 

Pconf 



8itc 2 Alal + 8 + 45p 2 b 2 a% + 16cA ^aJ - 24/2bcA ^4 - 24/dbal 

(22) 

and : 



p 4 {60irc 2 A 2 a 7 + 60a* + 47 25 fx 2 b 2 a™ + l20cA ^Tral° - 840cA Q pb^j na\ 7 - 840 pba 7 ] 

l£>Tic 2 Alal + 16 + 90/i 2 6 2 afj + 32cA)V 7ra o - 48 fibcA ^7val - 48//6ag 

(23) 

With the wavefunction (3)in equation (10) i.e. including both relavistic 
and confinement effect the charge radius p 2 con f +re i and convexity parame- 
ter C conf+rel are given by : 



p 2 a 2 (4 - 26) (3 - 26) 

P conf+rel ~ 4[X 2 ] 

and 

_ A| (6 - 26) (5 - 26) (4 - 26) (3 - 2e) [X 3 ] 

C conf+rel - gprj (25) 

where Xi,X 2 and X 3 are given in Appendix. 

We note that equations (25) and (26) are equivalent to equations (18) 
and (19)of Ref[19] exhibiting explicit cA dependence. 



3 Results 

3.1 Values of cAq and lower bounds on a s 

As noted earlier , cA Q depends on /i,Mp,fp and ct s .In fig.l(a-e) we plot cA 
vs a s for D,D S ,B,B S and B c mesons .It shows that a s tends to reach the 
minimum value when two solutions of Eq.(ll) almost merge satisfying the 
condition(16).This feature is true for any set of the parameters p, f p and M p . 
In table 1 ,we give the lower bounds on a s for mesons having c and b quarks. 

The dependence of cA on a s and /i can be noted as follows : 
With constant \i , cA decreases with a s values rising and vice- versa. On the 
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Table 1: Lower Bounds on a s 



Mesons 


Quark content 


At(GeV) 
Ref[25] 


M p (GeV) 
Ref[25] 


U(GeV) 
Ref[25] 


cA 


Lower bound 
on a s 


D 


cu/cd 


0.276 


1.869 


0.192 


0.9665 


~ 0.601 


B 


bu/bd 


0.315 


1.968 


0.157 


0.7653 


~ 0.652 


D s 


cs 


0.368 


5.279 


0.210 


1.1967 


~ 0.49 


B s 


bs 


0.44 


5.279 


0.171 


0.999 


~ 0.493 


B c 


bc 


1.18 


5.37 


0.36 


1.167 


~ 0.302 



other hand, with constant a s , cA increases(decreases) with increase (de- 
crease) in /i. 



3.2 Bounds on slope and curvature of the IW function 

Using the lower bounds on a s for each heavy-light and heavy-heavy mesons, we 
obtain upper bounds on the slope and curvature of the I-W function using 
equations (23), (24), (25) and (26). They are listed in table 2. We note that 
with increasing a s values, the slope and curvature decreases and henceforth 
the lower bound on a s corresponds to the upper bound on p 2 and C. 

In table 3, we record the predictions of the slope and curvature of the 
IW function in various models while in table 4, we reproduce the correspond- 
ing predictions of the model of Ref.[19] with c = lGeV and A Q = 1 in V- 
scheme [251 EZ] for various mesons. Two values for B, B s and B c mesons are 
shown where ) represents the actual values for p 2 and C in that work 

with a s = 0.261;while case-b) represents those for adhoc adjustable value of 
a s = 0.60 in order to show the usefulness of large a s as mentioned in Ref 
19. The a s are already large for D and D s mesons, so no two values are shown. 
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Figure 1: Variation of cAq vs a s for D Meson. The +ve (-ve) solution of Eq.15 
corresponds to the dashed (solid) line and the two lines nearly coincide at a s 
~ 0.601 , the lower bound on a s corresponding to the solution of Eq.16 for 
D Meson. 
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< 

o 



2.2 
2 
1.8 
1.6 
1.4 
1.2 



-ve solution 
+ve solution 



1 

0.8 
0.6 



0.4 



0.48 



0.485 



0.49 



0.495 



alpha , 



0.505 



0.51 



0.515 



Figure 2: Variation of cAq vs a s for D s Meson. The +ve (-ve) solution of Eq.15 
corresponds to the dashed (solid) line and the two lines nearly coincide at 
a s ~ 0.49 , the lower bound on a s corresponding to the solution of Eq.16 for 
D s Meson. 
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1.8 



1.6 



1.4 



1.2 



< 

o 



0.8 



0.6 



0.4 
0.652 



~i 1 

-vg solution 
+ve solution 



0.654 0.656 0.658 



0.66 0.662 0.664 0.666 0.668 0.67 
alpha s *" 



Figure 3: Variation of cAq vs a s for B Meson. The +ve (-ve) solution of Eq.15 
corresponds to the dashed (solid) line and the two lines nearly coincide at a s 
~ 0.652 , the lower bound on a s corresponding to the solution of Eq.16 for 
B Meson. 
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Figure 4: Variation of cAq vs a s for B s Meson. The +ve (-ve) solution of Eq.15 
corresponds to the dashed (solid) line and the two lines nearly coincide at a s 
~ 0.493 , the lower bound on a s corresponding to the solution of Eq.16 for 
B s Meson. 
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2.2 



1.8 



1.6 



1.4 



1.2 



0.8 



-ve solution 
+ve solution 



0.305 



0.31 



0.315 



alpha , 



0.32 



0.325 



Figure 5: Variation of cAq vs a s for B c Meson. The +ve (-ve) solution of Eq.15 
corresponds to the dashed (solid) line and the two lines nearly coincide at a s 
~ 0.302 , the lower bound on a s corresponding to the solution of Eq.16 for 



B c Meson. 
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Table 2: Upper Bounds on Slope and Curvature 



Meson 


Slope p 2 


Curvat 


ure C 


(Quark 


Without relati- 


With relati- 


Without relati- 


With relativi- 


Content) 


vistic effect 


vistic effect 


vistic effect 


vistic effect 


D(cu/cd) 


6.78 


1.675 


13.19 


5.138 


B{bu/bd) 


5.78 


1.016 


9.58 


1.29 


D s (cs) 


9.115 


3.067 


26.48 


14.32 


B s (bs) 


11.92 


2.652 


34.49 


6.902 


B c (bc) 


28.46 


10.39 


219.46 


45.23 



Table 3: Predictions of the slope and curvature of the IW function in various 
models. 



Model 


Value of p 2 


Value of curvature C 


Yaouanc et al [28J 


> 0.75 




Yaouanc et al [T2] 


> 0.75 


> 0.47 


Rosner et al [29] 


1.66 


2.76 


Mannel et al§0\ ED] 


0.98 


0.98 


Pole Ansatz [32] 


1.42 


2.71 


MIT Bag Model [H] 


2.35 


3.95 


Simple Quark Model [3] 


1 


1.11 


Skryme Model [15] 


1.3 


0.85 


QCD Sum Rule $T5\ 


0.65 


0.47 


Relativistic Three Quark Model [I] 


1.35 


1.75 


Infinite Momentum Frame Quark Model [5] 


3.04 


6.81 
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Table 4: Predictions of the slope and curvature of the I-W function in the 
QCD inspired quark model according to Ref[19] with c = 1 and A Q = 1 
taking relavistic and confinement effect in V-scheme.This table is nothing , 
but the copy of the last rows of tables 1,2,3 of Ref 19 . 



Meson 


a s 


slope (p 2 ) 


curvature (C) 


D 


0.625 


1.136 


5.377 


D s 


0.625 


1.083 


3.583 


B 


<z)0.261 
6)0.60 


a)128.128 
6)1.329 


a)5212 
6)7.2 


B s 


a)0.261 
6)0.60 


a)112.759 
6)1.257 


a)4841 
6)4.379 


B c 


o)0.261 
6)0.60 


a)44.479 
6)1.523 


a)2318 
6)0.432 



4 Conclusion and Remarks 

In this paper,we have shown that the reality bound on cA puts lower limit 
on a s and correspondingly upper limit on p 2 and C. 

Furthermore, with cAq , the upper bounds on p 2 and C decrease which is 
evident from the above list of bounds [tab le-2]. The estimated upper bounds 
on p 2 and C for all the mesons are found to be consistent with other mod- 
els and data [table-3] without making any adhoc enhancement of the strong 
coupling constant as had been done in ref (19) [table-4] .From the phenomeno- 
logical point of view we note that in the nonrelativistic limit ,the universal 
form factor and Isgur-Wise function for semileptonic decay B — > D*lv are 
identical when subleading terms in velocity and terms of order O (J^j are 
neglected with Eb as the binding energy and ttiq as the mass of heavy quark 
[33] .However even if we make calculation for the universal form factor for fi- 
nite mass, we obtain to first order in (y — 1) as 0.8 -2.57 (y — 1) which seems 
to be satisfactory [33| 134"]. 

It is worth notable that in the limit cA — > 0, there will be no bounds 
on a s as well as on p 2 and C ; rather fixed values of a s have to be used to 
get definite set of p 2 and CSo,in that case, the analysis will turn to that of 
ref [17,18] where large confinement could not be (i.e. 6 = 0.183GeV 2 ) incor- 
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porated e.g. tables -(1,3) of ref[17] and tables -(2,3) of ref[18]. 



We conclude this paper with a comment on the physical significance of 
the factor 'd that has become so crucial for our analysis of bounds on slope 
and curvature. 

It is common wisdom that a constant potential like V just scales the en- 
ergies and doesnot affect the wavefunction nor does it change physics. This 
can be seen from the hydrogen atom problem with the potential V(r) == 
— + c.However , if one uses V as the perturbation instead of as parent in 
the Dalgarno method of perturbation theory |20] .the wavefunction for the 
fT-atom becomes: 

ip (r) = N x (cA Q + -±=) e"^ 
V v """o / 

to be compared with the wavefunction with 'd as parent: 

ip(r) = ( — 7= J e~^o 
V V wa o / 

where the normalization constant: 



Nf 



Thus, the perturbative child 'd rather than the parent 'd plays the crucial 
role in the present analysis. 



References 

[1] M. Neubert,Phys Rep 245,259 (1994) 

[2] N Isgur and M B Wise Phys Lett B 232,113 (1989) 

[3] B.Holdom,M. Sutherland and J Mureika,Phys Rev D 49,2359 (1994) 

[4] Ivanov M A ,Lyubouvitskij V E,K6rner L G,Kroll P ,Phy Rev D 
56,348(1997) 



16 



[5] K6nig B,K6rner J G ,Kramer, Kroll P ,Phy Rev 56,4282(1997) 
[6] F E Close and A Wambach ,Nucl Phys B 412,169 (1994) 
[7] H W Huang,Phys Rev D 56,1579 (1979) 
[8] D Melikhov, Phys Rev D 53,2460 (1996) 

[9] M R AhmadyR R Mendel and J D Talman,Phys Rev D 52,254 (1995) 

[10] M G Olsson and S Veseli,Fermilab-Pub-96/418-T 

[11] M G Olsson and S Veseli,Phys Rev D 51,2224 (1995) 

[12] A Le Yaouanc ,L Oliver and J C Raynal ;Phy Rev D 69,094022 (2004) 

[13] Y B Dai,C S Huang,M Q Huang and C Liu;Phys Lett B,387,379 (1996) 

[14] M Sadzikowski and K Zalewski ; Z Phys C 59, 667 (1993) 

[15] E Jenkins ,A Manohar and M B Wise ;Nucl Phys B 396,38(1996) 

[16] D K Choudhury ,P Das ,D D Goswami and J N Sharma ; Pramana 
J.Phys.44,519(1995) 

[17] D K Choudhury and N S Bordoloi ;IJMPA,Vol.l5,No.23,3667(2000) 

[18] D K Choudhury and N S Bordoloi ; MPLA, Vol. 17,No.29, 1909 (2002) 

17 



[19] D K Choudhury and N S Bordoloi ;MPLA,Vol.24,No.6,443(2009) 

[20] A Ghatak and S Lokanathan in "Quanum Mechanics" ,McGraw 
Hill,1997,pp-291 

[21] E Eichten ; Phys.Rev.D 17,3090 (1978) 

[22] J J Sakurai in "Advanced Quantum Mechanics" 
(Massachusetts, Addison- Willey Publishing Company,1986) pp-128 

[23] C Itzykson and J Zuber in "Quantum Field Theory" (international 
Student Edition,McGraw Hill, Singapore, 1986) pp-79 

[24] V O Galkin,A Yu Mishurov and R N Faustov;Sov J Nucl Phys 
53,1026(1991) 

[25] D E Groom et al(Particle Data Group) ;Eur Phy J C,15, 1(2000) 

[26] Y Schroder ;Phys Lett B 447,321(1999) 

[27] Y Schroder ; Nucl Phys Proc Suppl 86,525(2000) 

[28] A Le Yaouanc ,L Oliver ,0 Pene and J C Raynal ;Phys.Lett.B 365; 319 
(1996) 

[29] J L Rosner ; Phys Rev D 44,3732 (1990) 

[30] T Mannel ,W Roberts and Z Ryzak;Phys Rev D 44,R18 (1991) 

18 



[31] T Mannel,W Roberts and Z Ryzak;Phys Lett B 255,593(1993) 



[32] M Neubert ;Phys Lett B 264,455 (1991) 



F Jugeau , A Le Yaouanc,L Oliver, J C Raynal;Phys Rev D 
70,114020(2004) 



[34] F E Close , A Wambach;RAL-94-041, OUTP-94 09P; April 1994 
Xi, X 2 and X 3 are evaluated as : 

X l = Q47rc 2 A 2 a 3 + 64 + /i 2 b 2 a 6 (8 - 2e) (7 - 2e) (6 - 2e) (5 - 2e) 



+128cA V7rag - 16cA //&V7ra§ (6 - 2e) (5 - 2e) 

-16fjba% (6 - 2e) (5 - 2e) (26) 



X 2 = 647rcMSaj; + 64 + //Vag (6 - 2e) (5 - 2e) (4 - 2e) (3 - 2e) 
+ 128cA v^rag - 16cAo^&0ra§ (4 - 2e) (3 - 2e) 

-16(j,ba 3 (4 - 2e) (3 - 2e) (27) 



X 3 = 64vrc 2 A^a^ + 64 + //Va° (10 - 2e) (9 - 2e) (8 - 2e) (7 - 2e) 



+ 128cA \/7ra 3 - 16cA 0/ u6^7ra|] (8 - 2e) (7 - 2e) 

-16^6ajj (8 - 2e) (7 - 2e) (28) 

Not only the above expressions ,but all the integrals in the analysis are eval- 
uated with the help of Gamma function given by : 

r (n + 1) [+™ 



a n+1 



r+oo 

/ r n e- ar dr (29) 
Jo 
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